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The insulating state of matter is characterized by the excitation spectrum, but also by qualitative 
features of the electronic ground state. The insulating ground wavefunction in fact: (i) sustains 
macroscopic polarization, and (ii) is localized. We give a sharp definition of the latter concept, 
and we show how the two basic features stem from essentially the same formalism. Our approach 
to localization is exemplified by means of a two-band Hubbard model in one dimension. In the 
nonintcracting limit the wavefunction localization is measured by the spread of the Wannier orbitals. 



In a milestone paper appeared in 1964 Jy W. Komi 
investigated the very basic features which discriminate 
between an insulator and a metal: he gave evidence that 
localization of the electronic ground wavefunction implies 
zero DC conductivity, and therefore characterizes the in- 
sulating state. In this Letter we provide a definition 
of localization which is deeply rooted into the modern 
theory of polarization (2j-[|], and rather different from 
Kohn's one. Indeed, besides zero DC conductivity, the 
property which obviously discriminates between insula- 
tors and metals is dielectric polarization: whenever the 
bulk symmetry is low enough, an insulator displays non- 
trivial static polarization. Here we show that the whole 
information needed for describing both localization and 
polarization is embedded into the same many-body ex- 
pectation value: namely, the complex number zpf de- 
fined in Eq. (|l^) below. It was previously shown j5j that 
macroscopic polarization is essentially the phase of Zn'- 
here we show that the modulus of z/v yields a definition of 
localization length which is sharper and more meaningful 
than the available ones. In our formalism a vanishing zjy 
implies a delocalized wavefunction and an ill-defined po- 
larization: this characterizes the metallic state. Our def- 
inition is first demonstrated for a one-dimensional crys- 
talline system of independent electrons, in which case 
our localization length coincides (for insulators) with the 
spread of the Wannier orbitals. We then study a two- 
band Hubbard model undergoing a Mott-like transition: 
both in the band regime (below the transition) and in the 
highly correlated regime (above the transition) the wave- 
function turns out to be localized, while the localization 
length diverges at the transition point, thus indicating a 
metallic ground state. Our approach to localization in 
a many-electron system sharply discriminates between a 
conducting and nonconducting ground state, yet avoid- 
ing any reference to the excitation spectrum. 

Let us start with a single one-dimensional electron: 
the distinction between localized (bound) and dclocal- 





FIG. 1. The distribution |t/j(:e)| 2 of a localized state within 
periodic Born-von-Karman boundary conditions 

ized (scattering) states is a clearcut one when the usual 
boundary conditions are adopted; much less so when 
periodic Born-von-Karman boundary conditions (BvK) 
are adopted, implying a ring topology for the one- 
dimensional system. Within the latter choice — which 
is almost mandatory in condensed matter physics — all 
states appear in a sense as "delocalized" since all wave- 
functions ip{x) are periodic over the BvK period: ip(x + 
L) = ip(x). We are going to show that the key parameter 
to study localization of an electronic state within BvK is 
the dimensionless complex number z, defined as 



dx e l ^ x \^(x)\ 2 



(1) 



whose modulus is no larger than 1. In the case of extreme 
delocalization one has IV^aOl 2 = V-^ an< ^ z = 0> while in 
the case of extreme localization 



\ip(x)\ 2 = S(x — xq — mL), 



(2) 



and we get z = e % ~t x °. In the most general case, depicted 
in Fig. [j], the electron density ^(x)! 2 can always be writ- 
ten as a superposition of a function n\ oc , normalized over 
(—oo, oo), and of its periodic replicas: 
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\ip(x)\ 2 = ^ ni oc (x - x - mL), 



(3) 



Both xq and n\ oc (x) have a large arbitrariness: we restrict 
it a little bit by imposing that xq is the center of the 
distribution, in the sense that /_ dx xn\ oc (x) = 0. 

Using Eq. (||), z can be expressed in terms of the 
Fourier transform of n\ oc as: 



i2l X n~ < 27T 

= e 4 ^ Xo ni oc {- — ) 



(4) 



The distinction between a localized and a delocalized 
state becomes clear if one studies the behavior of z when 
the BvK periodicity L is varied. For a localized state, 
in fact, the shape of n\ oc {x) is essentially L-independent 
(exponentially with L for large L), while the opposite is 
true for a delocalized state. If the electron is localized in 
a region of space much smaller than L, its Fourier trans- 
form is smooth over reciprocal distances of the order of 
L~ 1 and can be expanded as: 



n\ oc {-—) 



dxx 2 n loc (x) + <D(L~ 3 ). 

(5) 

Therefore at the increase of L, \z\ tends to 1 for a local- 
ized state, while it vanishes in the delocalized case. 

A very natural definition of the center of a localized 
periodic distribution |?/;(a;)| 2 is provided by the phase of 
z through the formula: 



(x) = — Im In z, 



(6) 



first proposed by Selloni et al. in Ref. Q to track the 
adiabatic time evolution of a single quantum particle in 
a disordered condensed system within BvK. The expec- 
tation value (x) is defined modulo L, as expected since 
^(a;)! 2 is periodic: the previous equations imply indeed 
(x) ~ xo mod L. The modulus of z can be used to mea- 
sure the localization length A. Using Eq. (||) we get 



In Id 



dx x 2 n\ oc (x), 



(7) 



and the spread of the electronic distribution can be de- 
fined through: 



lnlzl 



(8) 



which for large L goes to a constant limit for a localized 
state, and diverges for a delocalized one. Eq. (||) provides 
an alternative measure of localization with respect to the 
usual participation ratio Q . 



So much about the one-electron problem. We are now 
going to consider a finite density of electrons tlq: N par- 
ticles in a periodic box of size L. Eventually, the thermo- 
dynamic limit is taken: N — > oo, L — > oo, N/L = no con- 
stant. Even for a system of independent electrons, our 
approach takes a simple and compact form if a many- 
body formulation is adopted. In this case the ground 
state obeys BvK in each electronic variable separately: 



,x N ) = *(xi, . . . ,Xi+L, . 



,xn) 



(9) 



Spin variables are not explicitated (here and in the fol- 
lowing formulas), while of course are taken care of in 
the calculations. In analogy with the one-particle case, 
we define the many-body multiplicative operator X = 



52i=i x ii an d the complex number 



ZN 



(*|C 



(10) 



which will be used to discriminate between a localized 
many-body ground eigenstate (where \zff\ — > 1 for large 
N) and a delocalized one, where z^ vanishes. Ergo, fol- 
lowing Kohn's viewpoint |l|, the modulus of zn will be 
used here to discriminate between insulators and metals. 
We start with the dimensionless quantity 



D = - lim Nln\z N \ 

N — »oc 



(11) 



which is finite in insulators and divergent in metals: we 
define the localization length as A = \/D/(2irna). We 
emphasize that our definition of localization — as well as 
the definition of polarization given in Ref. [|— deals on 
the same ground with a general system, either ordered 
or disordered, either independent-electron or correlated. 

For a crystalline system of independent electrons the 
many-body wavefunction ^ can be written as a Slater 
determinant of Bloch orbitals and z^ factorizes. Using 
the same algebra as in Ref. ||, Eq. (14) onwards, one 
can easily prove that for a metal zn vanishes, while for 
an insulator D converges to the Brillouin-zone (BZ) in- 
tegral: 



D = 4m 6 — / dk [V (u' k \u' k ) 

JaL ra=l 
- { U 'm,k\ u l,k){ui,k\Um,k)}- 



(12) 



Lm—1 



In Eq. ([12]) we assume a linear system of lattice constant 
a with mb occupied bands and density no — 2mb/a; u m ^ 
is the periodic factor in the Bloch orbital (chosen to be 
a differentiable function of k), and the prime indicates 
the fc-derivative. The integral in Eq. ( |l2"| ) is a "geometric 
distance" || and measures the spread A 2 = (a; 2 ) — (a;) 2 
of the optimally localized Wannier orbitals |^,|l7J : our A 2 
coincides in fact with Qi/m,b, where fli is the quantity 
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FIG. 2. Solid lines: —N\n \zn\ 2 as a function of N at sev- 
eral values of A/t for independent electrons (U = 0). From 
bottom to top, A/t assumes the values: 5, 2, 1, 0.5714, 0.1, 
0.01, respectively. The dashed line is obtained replacing the 
logarithm with its leading expansion at A/t — 0.5714. 



defined by Marzari and Vanderbilt (in one dimen- 
sion). 

Next we study a one-dimensional model of a corre- 
lated polar crystal. We focus on the centrosymmetric 
case where zn is real and its phase is either or tt. We 
choose a two-band Hubbard model at half filling, whose 
Hamiltonian is: 

53 (-i) 3 1 A sW - t(c]*Cj+i* + H.c.) ] + U nn n n > 

i° i 

(13) 

and depends on two parameters besides the Hubbard U: 
the hopping t, and the difference in site energies 2A. 

In the special case U — we recover a system of inde- 
pendent electrons and the model describes an insulator 



whenever A^0. As discussed above, D is finite, Eq. (12) 



in the insulating case and formally infinite (even at finite 
N) in the metallic one. We show in Fig. || the convergence 
of D for several values of A/t: the localization length 
diverges upon approaching the metallic state (A = 0), 
and a large system size iV is needed for evaluating D if 
A/t is small. Approximating D with its finite- N value 
is exactly equivalent to a discretization of the BZ inte- 
gral in Eq. (|l2|): if we further replace the logarithm in 
Eq. ( |TT| ) with its leading expansion, we recover the same 
discretization proposed in Ref. |fl(i|| . While one gets the 
same limiting value, our logarithm form converges much 
faster: this is also shown in Fig. ^, for a selected value of 
A/t. 

The case of U > is much more interesting, since no 
Wannier functions or single-particle orbitals can be de- 
fined. Notwithstanding, our A mantains its value of a 
meaningful measure of the localization of the many-body 
wavefunction as a whole, even in the highly correlated 
regime. The model Hamiltonian of Eq. (fi~3) has been 
thoroughly studied by several authors |12|-|14[]: when U 
is increased to large values, the system undergoes an in- 
teresting transition, from a band insulator to a Mott insu- 




FIG. 3. Dimensionless localization parameter D, Eq. (ftl|), 
as a function of U/t, where A/t = 0.5714. Dashed line: cal- 
culations at N = 8, which are not converged since the size is 
too small. Solid line: values of D, where the thermodynamic 
limit is achieved by means of the ansatz wavefunctions (see 
text); the divergence at the Mott transition is perspicuous. 



lator. The real number Zjj changes sign at the transition: 
this fact has an important physical meaning, since it in- 
dicates a swapping of roles between anion and cation. In 
the low U regime the anion is the ion having the lowest 
on-site energy (odd j for positive A) , while the opposite 
is true in the highly correlated regime. We have shown 
in Ref. |i~3] | that the anion-cation swap manifests itself in 
a discontinuous change of the dynamical (or Born) ionic 
charge, while instead the static charge is continuous and 
carries no information about the transition. The many- 
body wavefunction is explicitly needed for detecting the 
transition, and the relevant information is indeed em- 
bedded in zjv. We adopt in the followingthe value of 
A/t = 0.5714, previously used in Ref. (jjfltl: the Mott 
transition occurs then at U = 2.27i. 

We perform exact diagonalizations for N=8 via the 
Lanczos algorithm, as described in Ref. JT3|. The re- 
sults are shown in Fig. |^, dashed line, and would indicate 
an increase of the localization length until the transition 
point, where a discontinuous drop occurs; in the high-{7 
region the wavefunction is strongly localized. However, 
upon performing the calculations in this way the finite 
size effects are clearly very relevant: this depends on the 
chosen value of A/t. Even at U=0 (where we can afford 
exact diagonalizations at arbitrarily large sizes) the value 
of D calculated at iV=8 differs from the fully converged 
value by 27% (see Fig. ||): in the correlated case the sit- 
uation is expected to worsen. We have performed a few 
calculations at different sizes, up to N = 12: the conver- 
gence turns out to be slow and oscillatory, with N = 4n 
and N = An + 2 following different trends. 

We overcome this drawback upon building approxi- 
mate wavefunctions for much larger sizes. At a fixed 
size N we perform several independent calculations, us- 
ing skew (quasiperiodic) boundary conditions with Bloch 
vector k over each electronic variable separately: 
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,Xn) 



IkL 



,X N ). 

(14) 



We choose M equally spaced values of k in the interval 
[0,2tt/L): 



2tt 

s ML 



,M-1. 



(15) 



Each of the <£'/ Cs is therefore BvK periodical over a pe- 
riod L' = ML, and we build an ansatz wavefunction 
for N' = MN electrons as the antisymmetrized prod- 
uct of the M different JV-particle wavefunctions $fe s . In 
the simple case of N = 1 this construction yields the 
Slater determinant of M orbitals, and is therefore the 
exact wavefunction for a system of N' = M nonintcr- 
acting electrons. Upon choosing N > 1 one allows the 
MN electrons to correlate, but only in clusters of N at 
a time: of course, our ansatz wavefunction has restricted 
variational freedom. At any given M, and for even N, 
the number z^i factorizes as [n5|: 



M— l 



zn' — zmn 



(16) 

For instance taking N = 4 and M = 3 the ansatz re- 
produces the exact 12 sites result (at U/t = 1) within 
3%. We then approximate the thermodynamic limit 
upon studying the large M limit at fixed N: 



D = - lim NM\n\z NM \ 



(17) 



The values of D calculated for M — 100 and N — 8 are 
plotted in Fig. |[ solid line: one clearly sees a divergence 
of the localization length at the Mott transition, while 
the wavefunction becomes localized again in the highly 
correlated regime. 

The many-dimensional generalization of the present 
formulation is not straightforward: its presentation is 
outside the scope of the present Letter. We mention here 
only a few main features, (i): A is essentially a unidimen- 
sional quantity, in the sense that one fixes a direction and 
defines a localization length in that given direction: say 
\ xx if we choose zn identical in form to Eq. (|Toj) . For 
anisotropic crystals, different A's coexist: for instance in 
graphite we expect A to be finite in the direction normal 
to the basal plane, and divergent in the planar direction, 
(ii): For a d-dimensional system of N electrons in a cu- 
bic box of volume L d the factor N appearing in Eq. (|ll| ) 
must be replaced with jV 2 / d_1 in order to define the di- 
mensionless D XXl and the localization length becomes: 
\ 2 xx = A 1 / d - 1 D xx /(2Trn 1 /d ) 2 . (iii): For a crystalline sys- 
tem of independent electrons X xx can be expressed as a 
BZ integral. In three dimensions, for a cubic lattice and 
mj) occupied bands, we can prove that: 



\ z = 

"xx 



4-2/3 



(2nn 



1/3- 



2 N 



lim N- 1/3 ln\z N \ 2 



3771b 



n h (is) 



where f?i is the BZ integral of Ref. pi[ . 

In conclusion, we have shown how to unambiguously 
measure localization in the ground state of a many- 
electron system. We have shown over a few examples 
how to discriminate between an insulator and a metal 
without actually looking at the excitation spectrum, sim- 
ply scrutinizing electron localization in the ground eigen- 
state. For the special case of an insulating crystal of 
nonintcracting electrons we measure nothing else than 
the localization of the Wannier functions, whereas in the 
correlated and/or disordered case our approach to local- 
ization is not related — to our best knowledge — to any 
previously known theory Jl6| . Our work opens the way 
to further advances and leaves several important issues 
open. We mention just a few of them: effects (possibly 
qualitative) of long-range interaction upon A; use of dif- 
ferent approximate many-body wavefunctions (such as 
e.g. quantum Monte Carlo); role of A in the context of 
Anderson localization in disordered systems. 
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